Theory of itinerant-electron ferromagnetism 
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A theory of Kondo lattices or a 1/d expansion theory, with d spatial dimensionality, is applied 
to studying itinerant-electron ferromagnetism. Two relevant multi-band models are examined: a 
band-edge model where the chemical potential is at one of band-edges, the top or bottom of bands, 
and a flat-band model where one of bands is almost flat or dispersionless and the chemical potential 
is at the flat band. In both the models, a novel ferromagnetic exchange interaction arises from the 
virtual exchange of pair excitations of quasiparticles; it has two novel properties such as its strength 
is in proportion to the effective Fermi energy of quasiparticles and its temperature dependence is 
responsible for the Curie- Weiss law. When the Hund coupling J is strong enough, the superexchange 
interaction, which arises from the virtual exchange of pair excitations of electrons across the Mott- 
Hubbard gap, is ferromagnetic. In particular, it is definitely ferromagnetic for any nonzero J > 
in the large limit of band multiplicity. Ferromagnetic instability occurs, when the sum of the two 
exchange interactions is ferromagnetic and it overcomes the quenching of magnetic moments by the 
Kondo effect or local quantum spin fluctuations and the suppression of magnetic instability by the 
mode-mode coupling among intersite spin fluctuations. 
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I. INTRODUCTION 

Itinerant-electron ferromagnetism has been a long 
standing and important issue since the advent of quan- 
tum mechanics. When the Hubbard model, for example, 
is considered in the Hartree-Fock approximation, a fer- 
romagnetic instability condition is given by 



Upo{0) > 1, 



(1.1) 



with U the on-site repulsion and po(0) the unrenormal- 
ized density of states at the chemical potential. Local cor- 
relations between electrons with parallel spins are accu- 
rately treated in any approximation as long as the Pauli 
principle is taken into account. In the Hartree-Fock ap- 
proximation, however, no local correlations between elec- 
trons with antiparallel spins are taken into account; para- 
magnetic states are disfavored so that magnetic states are 
relatively favored. The condition (LI) is never reliable. 
Local spin correlations should be accurately taken into 
account far beyond the Hartree-Fock approximation in 
any reliable theory, j-. 

In 1963, KanamoritI studied ferromagnetism in Ni in 
the T-matrix approximation, which is relevant for low 
densities of electrons or holes. He showed that a ferro- 
magnetic state is stabilized not only when Eq. ( |li] ) is 
satisfied but also when the density of states has a sharp 
peak within broad bands and the chemical potential is 
at the peak. Mielkp3, Mielke and TasakiB, Tasakilj, and 
Kusakabe and AokiO showed that the ground state of cer- 
tain multi-band models is ferromagnetic. In their models, 
one of bands is dispersionless or fiat and the density of 
states is similar to that of Kanamori's. These theories 
show that the existence of a flat band or flat bands is 
one of the most relevant conditions for itinerant-electron 
ferromagnetism. In many ferromagnets such as Ni, the 
density of electrons or holes is low and the chemical po- 
tential is at one of the band edges. One may argue that 
this condition must also be relevant. Call these two con- 
ditions flat-band and band-edge conditions in this paper. 

The ground state can be ferromagnetic even if neitheg 
of these conditions is satisfied. NagaokaO and ThoulessQ 



showed that when a single hole is introduced into the 
just half- filled system or the electron number is A/" — 1, 
with N the number of lattice sites, the ground state is 
a completely polarized ferromagnetic state in the limit 
of Upo{0) +00. Because spin waves of this state are 
abnormal and any extension of their arguments to low 
but nonzero densities of holes is difficult, it is doubtful if 
a ferromagnetic state is stabilized in the thermodynamic 
limit of small hole densities where the number of electrons 
is N{1 - S), with N +oo and S +0. pp 

In 1963, anothGs two distinguished papersElfl were pub- 
lished. HubbardB showed that a band splits into two 
bands, the lower and upper Hubbard bands. When 
the density of states of unrenormalized electrons is a 
Lorentzian shape of width F l/7rpo(0), for example, 
the local part of the Green function is given by 



^^G^(i£n,k) = 



1 



i£n^p-ea- E^(i£n) +^Fsgn(£n) 

1 - jn-a) 
iSn^p- ea + iFsgn(£n) 

i€n^ p- Ca-U ^irsgn{en) 



, (1.2) 



with T^cj{ien) the selfenergy, /i the chemical potential, 
Ca the band center, and (ria) the average number of 
electrons-|path spin a per site. When the Fermi-liquid 
relationHEll is made use of, Gutzwiller's resultB gives the 
mass enhancement factor of 



4>7 - 1/(1 - 



(1.3) 



with n = (n|) + (^|), so that the quasi-particle band- 
width is approximately given by l/^^po(0). Hubbard's 
result is relevant for single-particle excitations far from 
the chemical potential, while Gutzwiller's result is rele- 
vant for those in the vicinity of the chemical potential. A 
combined theory of Hubbard's and Gutzwiller's gives a 
three-peak structure of the density of states, Gutzwiller's 
band at the chemical potential between the lower and up- 
per Hubbard bands. 



The selfenergy obtained by Hubbard does not depend 
on wave-numbers k and the mass enhancement factor 
obtained by Gutzwiher neither depends on k. The two 
theories are within the single-site approximation (SSA). 

Treating the Hubbard or the periodic Anderson model 
in one of the best SSA's that include all the single-site 
terms is reduced to determining and solYing selfconsis- 
tently a single-impurity Anderson model.ES The Ander- 
son model is one of the simplest effective Hamiltonians 
for the Kondo probleni|:^The Kondo problem is almost 
completely understood;t2rt^ ^^f^ exact solution by 
the Bethe method was obtained.tll Then, many useful re- 
sults are available to clarify single-site properties in lat- 
tice systems. The most essential physics involved in the 
Kondo problem is the quenching of magnetic moments 
by local quantum spin fluctuatiops, so that the ground 
state is singlet and nonmagnetiolj although it is close 
to a magnetic quantum critical point. No magnetic in- 
stability occurs within the SSA, and it can occur only 
when intersite or multi-site effects are properly consid- 
ered. Multi-site effects can be perturbatively considered 
by starting from an unperturbed state constructed in the 
SSA. Because the quenching of magnetic moments by lo- 
cal quantum spin fluctuations is properly treated in such 
a perturbative theory, it must give a reliable magnetic in- 
stability condition. Such a perturbative theory is nothing 
but a theory of Kondo lattices. It is also nothing but a 
1/d expansion theory, with d being spatial dimensional- 
ity. Any SSA that includes all the single-site terms is 
rigorous for Landau's -Eermi-liquid states in infinite di- 
mensions {d +00). LJ Multi-site terms are of higher 
order in l/d except for Weiss' mean fields, which are re- 
sponsible for magnetic instability. 

The spin susceptibility of magnets obeys the Curie- 
Weiss law. Two mechanisms hp^e been proposed for that 
of itinerant-electron magnets:E3 the suppression of the 
spin susceptibility by the^piode-mode coupling among 
intersite spin fluctu^k>nsElrE3 and its enhancement by 
Weiss' mean fields .EllES When the density of states is 
almost constant in the vicinity of the chemical poten- 
tial, the suppression by the mode-mode coupling be- 
comes larger wittu ijacreasing temperatures; it gives the 
Curie- Weiss law.E2rE3 When the density of states has a 
sharp peak at the chemical potential, on the other hand, 
the suppression becomes larger with decreasing temper- 
atures; the mode-mode coupling plays a totally negative 
role in the Curie- Weiss law in the flat-band model.E^ 

Because the mode-mode coupling suppresses the spin 
susceptibility in any case, it cannot be a mechanism of 
magnetic instability. A scenario that the mechanism of 
the Curie- Weiss law is also responsible for magnetic insta- 
bility itself seems to be more reasonable than the scenario 
that the mechanism of the Curie- Weiss law suppresses 
magnetic instability. 

The mode- mode coupling is of higher order in l/d so 
that it vanishes in the large limit of d. There must be an- 
other mechanism, which is of leading order in 1/d. Only 
the possiblej-ope is the temperature dependence of Weiss' 
mean fields.E2l The enhancement of the spin susceptibility 
by an exchange interaction arising from the virtual ex- 
change of pair excitations of quasi-particles can give the 
Curie- Weiss lawEj Polarizations multiplied by the ex- 
change interaction are Weiss' mean fields, and they are 
of leading order in 1/d. When the flat-band or band- 



edge conditions is satisfied, this exchange interaction is 
ferromagnetic and the Curie- Weiss law appears only in 
the homogeneous susceptibility. Miyai and the present 
authoiEj examined the Curie- Weiss law of Ni by using 
Kanamori's model. 

Field-induced ferromagnetism or metamagnetism is 
observed in many metals. In CeRu2Si2, for example, 
such a ferromagnetic state is stabilized at magnetic fields 
higher than Hm — 7.7 T. One of the most cruciaJLexperi- 
mental results is the single-parameter scaling ;E£lEZl many 
quantities sensitively depend on pressures, and they scale 
with the Kondo temperature or energy, ksTK- The 
Kondo energy is approximately equal to a half of the 
bandwidth of quasiparticles, ksTx — l/[2^^po(0)]. This 
experiment gives a strong restriction for a ferromagnetic 
exchange interaction responsible for the metamagnetism 
in CeRu2Si2. Satoh and the present author showed that 
the exchange interaction responsible for the Cuifie- Weiss 
law is also responsible for the metamagnetism.E^I In 4/ 
electron systems, the density of states is of a camelback 
structure or of a two-peak structure because of the hy- 
bridization between dispersive conduction bands and al- 
most dispersionless / bands. In the absence of magnetic 
fields, the chemical potential is between the two peaks; 
either of the flat-band and band-edge conditions is not 
satisfied, so that the exchange interaction is antiferro- 
magnetic. In the presence of fields as high as Hm-, the 
chemical potential is at one of the peaks of the camelback 
because of the Zeeman shift of up and down spin bands; 
both of the flat-band and band-edge conditions are sat- 
isfied, so that the exchange interaction becomes ferro- 
magnetic around Hm • Because this exchange interaction 
has a novel property that its strength is proportional to 
the bandwidth of quasiparticles l/(^^p(0), it is easy to 
explain the observed single-parameter scaling. Call this 
one simply a novel exchange interaction in this paper. 

The previous papersSj'ESi imply that when either or 
both of the flat-band and band-edge conditions are sat- 
isfied ferromagnetic instability is caused by the novel ex- 
change interaction even in the absence of fields. Howe^ie^, 
there is a crucial drawback in the previous paper s.SES 
An additional phenomenological ferromagnetic exchange 
interaction, which is independent of temperatures and 
magnetic fields, was assumed to obtain a quantitative 
agreement between experiment and theory. The authors 
of these papers speculated that such a ferromagnetic ex- 
change interaction must exist in multi-band models. 

The purpose of this paper is to study ferromagnetism 
in multi-band models from a theoretical viewpoint of the 
competition between magnetic instability caused by in- 
tersite exchange interactions and the quenching of mag- 
netic moments by the Kondo effect. This paper is orga- 
nized as follows: In Sec. |l|, an unperturbed state in the 
1/d expansion is constructed in the SSA. It is shown in 
Sec. |l|that when the Hund coupling is strong enough the 
superexchange interaction, which arises from the virtual 
exchange of pair excitations between the lower and up- 
per Hubbard bands, becomes ferromagnetic. It is shown 
in Sec. IV that when the flat-band or band-edge condi- 
tion is satisfied a ferromagnetic state can be stabilized. 
Discussion is given in Sec. and summary is given in 



Sec. VI 



In Appendix a single-impurity Anderson 
model is studied by the 1/(2/ + 1) expansion method, 
with 2/ + 1 being orbital degeneracy. 



II. PRELIMINARIES 



with 



A. Multi-band model with the Hund coupling 

Consider ^-.(2/ + l)-band model, with / an integer or a 
half integer :El 



1~L = TYo + 

The first term is a one-body term given by 



(2.1) 



^0 = [^£^mm^ (k) — jlSmm' ^^Jnkcr^^'ka , (2.2) 

mm' ka 



with a]^kcr ^ creation operator of carriers, electrons or 
holes. Bands or orbitals are denoted by m = — —(/ — I), 
• • •, and and up and down spins by cr = ±1. It is 
assumed that all the bands are degenerate at the T point 
due to lattice symmetry, so that both £^^^(0) and 



a — ^ Emmi}^) 



(2.3) 



do not depend on m, and E^^f (0) = and (l/VN) Ek 
Emm'0^)XinvO^) = for m 7^ and any function 
XinvO^) invariant under any symmetry transformation. 
In this paper, Emm'0^) = Em'mO^) being real and 



Err 



'(k) = 



£;;;(k), m 



(2.4) 



£^;_((k), m^m' 

are also assumed.0 This assumption makes arguments in 
this paper quite simple, as is shown later. The dispersion 
relations of 2? + 1 bands are given by 

_ ^ r Eu{k) + 2lEi-i{k), m = l 
E^{k) = <^ ' . 2.5 

\Eu{k)-Ei.i{k), -l<m<{l-l) 

The fiat-band condition is satisfied when the Ith band is 
almost dispersionless in such a way that 



Et(k) = Eu{k) + 2lEi_i(k) Eu{0) 



(2.6) 



for any k. In this case, the dispersion relations of other 
21 bands {m ^ I) are given by 

E*^{k.) Euik) + i [Eaik) - EuiO)]. (2.7) 

The second term in Eq. ( |2.lD is the on-site Coulomb 
interaction: 

mm' iaa' 

+ 1"^ Y.^^^iym'ia''^mia'am'i<T (2.8a) 

m^ 



m^m' iaa' 

^mia^mi— a 



1 / 1 \ 

'2 ( ^ 2^/ ^ ^ ^ '^mia'^m'ia' 
m^m' icTcr' 

m^m' iX a-it72 a-3<T4 

^ '^Licri <^mio-2 '^In'ias '^m.'iui ■, 



(2.8b) 



i(k-R.)^t 



^m\^(j^ 



(2.9) 



amia its Hermite conjugate, and Umia = ^^Lia^^^^- 
cause of the degeneracy of 2/ + 1 bands at the T point, 

= <^ , , (2.10) 

is also assumed. In general, the H und coupling is fer- 
romagnetic (J > 0). In Eq. ( p.8bD , 5^^^' is the (aa^th 
element of the ^S* = ^ spin matrix, Sy, or s^. In this 
paper, is diagonalized so that = ^aS^jcr' - 

B. Mapping to the Anderson model 

The single-particle Green function for H is given by 

g-\ien, k) = {iSn + /i) X - S{k) - S^iiSn, k), (2.11) 



with X being the unit matrix, [f(k)]^ 



= Er, 



'(k), 



and 5cr(^£n,k) the selfenergy matrix. The selfenergy is 
divided into single-site and multisite terms: 



(2.12) 



with T^(j{ien) the single-site selfenergy and ASa{iSn^'k) 
the multisite one. In our model, the selfenergy for k = 
is proportional to the unit matrix because of the assumed 
symmetry; [AScr{iSn^O)]^^ does not depend on m and 
[ASa{i€n,0)]^^, = for m 7^ m\ 

Calculating the single-site selfenergy T^criiSn) is re- 
duced to solving a single-impurity Anderson model: 



Ha 



[Ca - jj] rima + ^ [ec(k) - jj] cj^^Ck, 
ma kcr 

f r t * t 

/at / ^ r^mk^^rncr^kcr + '^mk^kcr^^ 



T,kcr 
mm' aa' 



Umm' (^rna^m' a' ^mcr 



Y1 "^^Ina^^ii'a'^rna'am'a. (2.13) 



m^m' crcr' 



with rima = cilna^mcT' Here, ea, /i, t/mm', and J are 
the same as those for the lattice system; ec(k) and Vm\<i 
should be determined so as to satisfy a mapping condition 
discussed below. The Green function for Ha is given by 



(2.14) 



with L^iiSn) = (l/A^)Ekbmk|V[^^n - ec(k)]. Be- 
cause of the assumed symmetry, Q~^{i£n) is diagonal and 
L(j{^£n) does not depend on m. Because the symmetri- 
cal property of [AtScr(^6:n, k)] is the same as that of 
^rnm'(k), the local Green function defined by 



^a{iSn) = ^ ^ga(^gn,k) 
k 



(2.15) 



is also diagonal. When 



(2.16) 



is satisfied any single-site quantity of the lattice system 
is the same as its corresponding quantity of the single- 
impurity Anderson model. The quantities, ec(k) and 
^mk, of the Anderson mod el sh ould be selfconsistently 
determined to satisfy Eq. ( 2.16| ); Eq. ( 2.16 ) is nothing 
but the mapping condition to the Anderson model. 



III. EXCHANGE INTERACTIONS 



A. Kondo lattices 



Consider the mapped Anderson model ( 2.13| ), and in- 
troduce fictitious and infinitesimally small external fields 
such as Hext = - Y^rna A^^^n^^, with 



A/i + ijlb 



^cr9sHs^mgoHo^^{crm)gsoHs, 



(3.1) 



Here, is the Bohr magneton; ^g, go and gso are g 
factorsEil The single-site selfenergy at T = K is ex- 
panded in such a way that 

^ma{i£n) = ^(0) + (1 - 4)^)iSn + (1 - ^c)A/i 

+ (1 - 4)s)^agsfiBHs + (1 - ^i)mgoliBHo 



+ (1 - (l)so)-{cFm)gsol~iBHs. 



(3.2) 



for \en\ <C ksTx^ with Tk the Kondo temperature dis- 
cussed in Introduction and defined below. 

When polarizations of conduction electrons are ig- 
nored, it follows according to Yoshimori'sE3 that 



m^^)k 

2(21 + 1) ^ 2(2/ + 3) (2/ + 1)(2/ + 3) 



, (3.3) 



and static susceptibilities at T = K of the Anderson 
model are simply given by 



Xc(0) 



d 



(3.4) 



o(0) -Q^YrngoliB{nma) 

^ ma 



(3.6) 



and 



XsoiO) 



dHs 



"X! \ i'^^)9sotJ'B{nm.a) 



(3.7) 



with 



x:(o) = 2(2; + i)0„/,(o), 



(3.8) 



for charge (a = c), spin (a = s), band or orbital (a = o), 
and spin-orbital combined (a = so) channels. Here, (• • •) 
stands for the thermal average, and 



P{e) = 



1 



2^ + 1 



-- )Im 

TT 



Tr g^ie + iO) 



(3.9) 



is the density of states per spin and band. 

So far, there is no restriction for [/, U' ^ or J. When 
U = U' and J = 0, = = ^so. It is reasonable 
to assume that U ^ J ^ U' ^ J and U ^ U' so that 
Only a single energy scale, the Kondo 
temperature, appears in low-energy phenomena of the 
single-site problem. In this paper, it is defined by 



21 -\- 1 



(3.10) 



The specific heat at low temperatures, T <^ Tk, is given 
by (7 = 7T + • • •, with 



(3.11) 



The Kondo temperature, ksTK, is as large as the effec- 
tive Fermi energy of quasiparticles in the lattice system. 

In this paper, the Kondo temperature is treated as a 
phenomenological parameter instead of determining and 
solving the Anderson model selfconsistently. 

Define a two-point polarization function of the lattice 
system, which is a 2(2/ + 1) x 2(2/ + 1) matrix, by 

\m^u q)] , , = - y / dr e-^-+^^-(i^^-i^^) 

L V ''^ \m(j,m'(T' i\T / ^ 

ij 

x{nmia{r)nm'ja'{0)). (3.12) 
When its irreducible part is denoted by 7r(zcc;/,q), 



n(icj/, q) = 7r{iu;i, q) [1 + UiriiuJi, q) 



(3.13) 



According to Eq. (|2.8b|) , U is defined by 



[U] 



ma,m' a' 



0, m = m' and cf = a' 

U' — m' and cj = a' 

[/, m = m' and cj ^ cj' 

U' ^ m ^ m' and a ^ a' 



(3.14) 



The irreducible part is divided into a single-site term 
T^{iuJi) and a multisite term A7r(za;^,q) so that 



7r(zcj^, q) = T^{iijJi) + A7r{iuJi, q). 



(3.15) 



The reducible polarization function of the mapped An- 
derson model is given by 

U{iiJi) = [1 + Un{iui)] ~\ (3.16) 

Define the following 2(2/ + 1) x 2(2/ + 1) matrix: 



(a) 



ma.m'a' 



2(2/ + 1) 
1 



(7(7 , a = s 



(3.17) 



a = o 



2(2/ + 1) 
1 3mm' 
2(2/ + l)/(/ + l) ' 
1 Smm'crcr' 
I 2(2/ + 1) /(/ + 1) ' 

Because our syst em has the particular symmetry such as 
Eqs. (|^, (|2.1Q|) and (|3.14D , susceptibihties of the four 
channels are simply given by 

xliiu^u^) ^ 2{2i + i)EE^itn.v [n(i-'<-q)]^.,^v 

ma m'a' 

TTaiiuJi.q) 



and 



with 



2(2/ + 1) 



1 - Ua7Ta{iu;i,q) 



(3.18) 



(a) 



tl{iuJi) 



ma,m' a' 



ma m'a' 



2(2^ + l)-4M_^, 

1 - UaTTaitOJl) 



(3.19) 



instead of Eq. (LI). When either of the flat- band and 
band-edge cond itions is satisfled, Eq. ( 3.26| ) is satis- 
fled. Eq. d^) leads to a set of UaXU^) > 1 and 
UaXai^^^) ^ 1 except for a = c. Then, as long as 

\cOl\<kBTK, 



Uana{iu;i) = 1^0{l/UaXa{^)) 

and Ua7Ta{iuJi,q) = 1 0{l/UaXa{^^^)) so that 
2 



(3.27) 



(3.28) 



The susceptibility given by Eq. (|3.24|) is consistent with 
a physical picture for Kondo lattices, where local fluctu- 
ations at different sites interact with each other by in- 
tersite exchan ge int erac tions . Then, we call Ia{iuJi^q) 
given by Eq. ( 3.25 ) or ( 3.28 ) an intersite exchange in- 
teraction for the spin, orbital, or combined channel. A 
perturbative treatment of Ia{iuJi^q) is nothing but the 
1/d expansion; Ia{iooi^ q) is of higher order in 1/d for al- 
most all q's, and it is leading order in 1/d for particular 
q's such as q = 0, nesting wave vectors, and so on. 



/ J / J ma.m'a' \_^^\raa,m'a' 
ma m' a' 

( -U -2l(2U' - J), a = c 
U + 2/ J, a = s 

-U^2U'-J, a = o 
U — a = so 



(3.20) 



Here, 7ra{i^h(l) and na{i^i) are irreducible polarization 
functions deflned by 



^aiiu^l) =^Y1 ^mlm'a' ^(^^0 



ma m'a' 



J ma,m' a' 



(3.21) 



A7ra(iw;,q) =^^ii:^j_„,^, [^A7r(iw,, q; 

ma m'a' 

and 

TTc {iuJi , q) = 7fc (zcJz ) + AtTc {iuJi , q) . 



ma,m'a' 



, (3.22) 



(3.23) 



Note that Xa{iu;i) is nothing but the susceptibility of the 

mapped Anderson modeL 

It follows from Eqs. (|3l8| ) and i^J^) that 



Xa(^^/,q) 



1 - |/c,(zcj^,q)x*(icj^)' 



with 



Ia{iu;i,q) 



ATTaiiuJl.Cl) 



(3.24) 



(3.25) 



2/ + 1 7Ca{iu;i)7ra{iu;i,q) 
In this paper, the strong coupling regime is deflned by 
Us > ksTK (3.26) 



FIG. 1. Examples of two-line diagrams for Att^ (zu;^ , q) . 
Diagram (a), which is in the site representation, gives the 
superexchange interaction between the ith and jth sites. Di- 
agram (6), which is in the wave- number representation, gives 
the novel exchange interaction, which arises from the virtual 
exchange of pair excitations of quasi-particles with wave num- 
bers k + q and k. Two-line diagrams higher order in intersite 
effects such as (c) can be still of leading order in 1/d, but 
they are ignored in this paper. A solid line stands for the sin- 
gle-particle Green function, Gma{i£n) given by Eq. (3.31) in 
(a) or the intersi te pa rt of G(r{i£n, k+q) or Ga(iSnj k) given by 
Eqs. ( |3.49| ) and (Oo|) in (b). A cross in (a) stands for tmi,m'j 
given by Eq. (|3.38D, a shaded triangle for the single-site ir- 



reducible t hree- point vertex function given by Eq. (3.34) in 
(a) or Eq. (3.51) in (6), and a shaded rectangle in (c) for the 
single-site irreducible four-point vertex function. 

In general, an exchange interaction arises from the vir- 
tual exchange of bosons or boso nic e xcitations. The ex- 
change interaction given by Eq. ( 3.28 ) is divided into four 
terms in such a way that 



la [itoi , q) = {iiOi , q) + 1^^^ (iw; , q) 

+Ali''^HiiOi,q) - 4A„(ic^,,q). (3.29) 

Only contributions from the so called two-line diagrams 
can be of leading order inl/d. Three examples of two-line 
diagrams are shown in Fig. 0. The first term, /^^^ {icui, q), 
corresponds to Fig. ||(a), and is due to the virtual ex- 
change of pair excitations of electrons or holes across 
the Mott-Hubbard gap, whose excitation energies are as 
large as or U' ± J and much larger than ksTK- This 
is a generalized superexchange interaction. The second 
term, I^\iui^c\)^ corresponds to Fig. and is due 

to the virtual exchange of pair excitations of quasipar- 
ticles, whose excitation energies are ©(/c^Tx). This is 
the novel exchange interaction discussed in Introduction. 
The third term, A/^^^^ (zcc;^, q), is the contribution from 
two-line diagrams higher order in intersite effects such as 
Fig. 0(c). The last term, — 4AcK(za;^, q), includes all the 
other contributions, which are definitely higher order in 
1/d. Its lowest-order Feynman diagrams correspond to 
those of the mode-mode coupling terms considered in the 
so called selfconsisteat renormalization (SCR) theory of 
spin fluctuations .E3'E3 

The super excha nge interaction, /^^^(zcj^, q), is consid- 
ered in Sec. Ill B , and the novel exchange interaction, 
I^\iuJi^q)^ is considered in Sec. IIIC. On the other 
hand, Aa{iuJi^c[) as well as AIa^\iui^ci) are ignored in 
this paper. Because AQ,(icj^,q) is of higher order in l/d 
for any q, it can be ignored at least for large enough d. 



B. Superexchange interaction 

There is no restriction for carrier densities so far. In 
the following part of this jpaper, our study is restricted 
to low densities of carriers^ 



{rima) < 1, ri = y^^jrima) < 1, 



(3.30) 



so that multiply occupied configurations at a single site 
can be ignored in the ground state and low lying excited 
states. The single-particle Green function is approxi- 
mately given by [Ga{i£n)]^^, = ^mm'Gma{i£n), with 



^ (j^m — cr) 

i£n^ l-i- ea-U ^ irsgn(£n) 

{rim' a) H~ 2{'^m' — a) 



-E 



2 (P"m' — a) 



i£n^ IJi- ea-U' - J ^ iVsgli^Sn) 



Here, ttF :^ l/po(0) and 



1 - ^{rim'a') 



(3.31) 



(3.32) 



is the density of empty sites. This is a stirjiightforward 
extension of Eq. (|L2|) or Hubbardis. theoryfl 

According to the Ward identity,LJ the reducible single- 
site three-point vertex function is given by 



J ma,m a 



dAEm'a' ' 



(3.33) 



When only t he de pendence of {rima) on AEm'a' is consid- 
ered in Eq. ( 3.31 ), the irreducible single-site three-point 
vertex function is calculated in such a way that 

ma m'a' 

X [A(i£„,ie„;0)[l-iYfi-(0)] 



ma,m a 



E 



An 



7ra(0) 

G^^(i£„) ^ isn -Zr + irsgn(e„) ' 

(3.34) 



with 



( - ^ ^^ 

ea^U' -J -II 



(3.35) 



and 



( As,i Ao,i Aso:i\ I 1 

^s;2 Ao-2 Aso:2 _ ~1 

As-3 ^o;3 ^so:3 I 

\ ^s;4 Aq-4 Asq-4 j 



1\ 

-1 

1 
" 2 



(3.36) 



for spin (a = 5), orbital (a = o), and combined (a = so) 
channels. 

Consider an exchange interaction between the ith and 
jth sites, which is defined by 

= \Y.^^\i^u^)e-'^-'^^^-^^ (3.37) 
q 

with being the position of the ith unit cell. Trans- 
fer integrals between_the ith and jth sites have also the 
following symmetry 



t 



m = m 



(3.38) 



When the two-line diagram shown in Fig. |i|(a) is taken 
into account and Eq. ( |3.27[) is made use of, the superex- 
change interaction is calculated in such a way that 



2/ + 1 U 



TT a 



(3.39) 



with 



/ J -^^ma,m'a\^T^'^^ 



m'j I 



' 1 2 



P, a = o and so 



(3.40) 



and 

ya;ij = \u{-kBT)^Uf\l{ien,ien;0)Gi,{ien) 

en 

1 ^ 

= 2 ^ ^ ^a;rA„;^/S(Zr, Zr'] T). (3.41) 

Here, z'; F) is defined in such a way that 



s(z,z^r) = 



Tl{z — Z') 



Tan"^ 



-Tan"^ 



(3.42) 



for z 7^ z', and = hm^/^^ z''). Note that 

= — 1 for r/[/ ^ +0 a nd J = 0. When the carrier 
density is so smah that Eq. ( 3.30 ) is satisfied, the contri- 
bution from excitations within the lower Hubbard band 
(r = = 1) is nothing but that from the virtual exchange 
of pair excitations of quasi-particles. Because it is con- 
sidered in the next subsection, it should be subtracted 
in order to -avoid double counting so that Eq. (3.41) is 
replaced byEZi 



(3.43) 



In the following part, I^.lj{iiJi = 0) is simply denoted by 

I^-lj. Because its uji dependence is small, it can be used 
for \uji\ <C/7and - J. 

In the small limit ofT/U +0, it follows that 



As) 

s;ij 



2/ + 1 



1 

U 



1 



1 



U'-J U'^J 



(3.44a) 



the superexchange interaction is ferromagnetic even for 
non-zero F. 



FIG. 2. Ys-ij defined by Eq. (3.43) as a function of 
the Hund coupling, J, for / = 1 and / = 2. Here, 
{ta — ii)/U — —0.1 and U' /U = 3/4 are assumed. A solid line 
shows Ys-ij for /U = 0.3, a broken line for nV/U = 0.2, 
a dotted-broken line for ttT/U = 0.1, and a dotted line for 
nV/U = 0. 
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4T, 



2/ + 1 
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U 



(3.44b) 



and 



As) 

so:ij 



4T.^ 



(so) 



2/ + 1 



1 

U 



U'-J U'^J 



(3.44c) 



In deriving these results, Ca — jii <^ -ksT^ Ca ^U — ji^ 
ksT and Ca ^ U' — J — (i ^ ksT are assumed. These 
results agree with those calculated for in^ating phases 
by the conventio nal pe rturbative method.L^ 

Note that Eq. ( ^.44a ) can be applied even to / = 0. The 
superexchange interaction in spin channels can never be 
ferromagnetic in the single-band model (/ = 0). It can 
be ferromagnetic in multi-band models (/ > 1/2), when 
J is as large as 



J> \^VU^ ^iU') -lU. 



(3.45) 



It is interesting that as long as J > the superexchange 
interaction is definitely ferromagnetic in the large limit 
of orbital multiplicity (/ +oo). 

Non-zero F makes the superexchange interaction re- 
duced. Fig. H shows Ys-^ij as a function of J for various 
F's, and Fig. || shows Ys-ij as a function of ea — /i for 
the same F's. Positive and negative Yg-ij''^ mean that the 
superexchange interaction is ferromagnetic and antiferro- 
magnetic, respectively. When J and / are large enough. 



FIG. 3. 

depth, Ca 



■ij defined by Eq. ( p.43| ) as a function of the 
, for / = 1 and / = 2. Here, J/U = 0.3 
and U' /U = 3/4 are assumed. A solid line shows Ys-ij for 
ttT/U = 0.3, a broken line for ttT/U = 0.2, a dotted-broken 
line for 7rF/[/ = 0.1, and a dotted line for ttF/U = 0. 



The signs of I^^^j and Ig^^^j depend on \t 



is) 



When It 



so:zj 
is) 



It;.. 



both of ^7, and are small. In- 



teresting is the case where either of \tij\'^ and 
much larger than the other, because it is possible that an 



instability of the orbital or combined channel occurs at 
low enough channels. In general, \Io% \ ^ l^io-ijl- 



C. Novel exchange interaction 



N 

En k ram' a 

(3.53) 



Selfenergy corrections from multi-site effects are of 
higher order ml/d except for those due to Weiss' mean 
fields, which are not present in Landau's normal Fermi- 
liquid phases. When only the single-site selfenergy given 
by Eq. (3^) is taken into account, the single-particle 
Green function is approximately given by 



[Qa ^ (^^n, k)] = [iSn^mm' - ^mm' (k)] , (3.46) 



with 



ew(k) = J- [^mm'(k) + fEo(O) - /i) ^wj , (3.47) 

for \en\ ^ ksTK- The dispersion relation of quasiparti- 
cles, ^m(k), is determined from det [^"-^(^^(k), k)] = 0, 
so that 



^m(k) 



l]0(0)+^:,(k)-/i 



(3.48) 



with ^^(k) given by Eq. (^. Note that ^,^(k) = ^-^k) 
for m such as — / < m < (/ — 1) and only ^^(k) can 
be different from them. The (mm' )th component of the 
Green function, is calculated in such a way that 



21 



1 



(2/ + !)^^ \ien- ^-Kk) i^n- 6(k). 

(3.49) 

for intraband components, and 



(-1) 



21+1 



1 



1 



(2/ + l)(^^ [isn- ^-z(k) isn- 6(k). 



(3.50) 



for interband components (m ^ tti^. 
According to the Ward identity)^ 



A«(0,0;0) = [1- /7«^«(0)]^« 
_ 2(2/ + !)^, 



1 + 



(3.51) 



When this is approximately used for non-zero energies 
such as l^nl ^ ksTk and \uoi\ ^ ksTk and the two-line 
diagram shown in Fig. |l|(6) is taken into account, the 
exchange interaction due to the virtual exchange of pair 
excitations of quasiparticles is calculated so that 



j(Q)(iu;,,q) = 



C/„Aa(0,0;0) ^i- P«(iwi,q)-P«(iwi) 



8(2/ + 1) 



(3.52) 



Here, 



is a a polarization function of quasi-particles and 

PaiiuJi) = ^^Pa{iooi,q) (3.54) 



is its local part. Because any single-site term is con- 
sidered in the un pertur bed state, the single-site part is 
subtracted in Eq. ( 3.52| ). 

For the spin channel {a = 5), Eqs. ( 3.53 ) and ( |3.54 ) 
become so simple as 



Ps{iu;i,q) = 



/[U(k+q)]-/[^^(k)] 
(2/ + l)7Vf^ e^(k)-^^(k+q)-zc^, 



(3.55) 



and 



00 pOO 



Ps{0) = I dxl dy^^^^^^^p*{x)p*{y), (3.56) 



CO ^ — CO 



x-y 



with 



P*(e) = ^E'^(^-^-(^)) (3.57) 



the density of quasi-particle states per spin and band. 
For orbital {a = o) and combined {a = so) channels, 

P U , r.^ - ^ /fe(k+q)]-/fc(k)] 

81 ^ /[e-Kk+q)]-/[6(k)] 
{2l + l)N^ 6(k)-e-i(k+q)-iw( 



2^(2^-l) ^ /[e-Kk+q)]-/K-/(k)] 
^(2^ + l)Ar+ e-/(k)-e-/(k+q)-ia;r ' 

This exchange interaction includes no intrasite part 
so that Eq^«'^^(^^^q) = 0- When lP{0,^0) < 0, 
^Q;^^(0,q) takes its maximum value at non-zero q. In 
this paper, it is classified as being antiferromagnetic 
when (0, +0) < and as being ferromagnetic when 
I^\o^ +0) > 0. As was studied in the previous paper ,0 

the temperature dependence of /i^^(0,q) is responsible 
for the Curie- Weiss law. It is easy to see from Eq. ( ^.52D 

or (O) that li^\0, q) is of the order of ksTK. 



IV. FLAT-BAND AND BAND-EDGE MODELS 

A. Instability condition 

Assume T = K in this section. Then, it follows that 
^.(0,|q| ^+0) = p*(0) and 



i/i^)(o,+o) = t 



p*(0) 



1 



(4.1) 



Note that p*{0) = 4>^p{Q) to leading order in 1/d. An 
instability condition against ferromagnetism is simply 
given by 1 - (1/4)1,(0, 0)xs(O) < or 



p*{0) 



AX, < 0, 



with 



AX, 



/W(0,0)-4A, (0,0)1 x.(0). 



(4.2) 



(4.3) 



When J is strong enough, the superexchange interaction 
/i^''(0,0), is ferromagnetic, as is discussed in Sec. pIB 



The mode- mode couphng, A5(0,0), can be ignored in in- 
finite dimensions {d ^ -\-oo). It is hkely th at A Xg > 
when / and J are large enough. When Eq. (O) is posi- 
tively large enough, ferromagnetic instability occurs even 
when neither of the fiat-banJ-and band-edge conditions 
is satisfied. The SCR theoryElrE3 can only be applied to 
such ferromagnets in finite dimensions {d <C +00). Then, 
it is examined in this section whether or not ferromag- 
netic instability is possible even when AXg is as small as 
I AX, I < ksTK or AXs = 0. 

The signs of Eq. or [l - Ps{0)/ p'^iO)] are func- 

tionals of p*(e). In order to examine the fiat-band and 
band-edge conditions, we consider the following model of 
y9*(e): It is given by the sum of a narrow triangular one 
and a wide square one and the chemical potential is just 
at the peak of the triangle in such a way that 



p*(e)=cplie) + il-c)plie), 
where p* (e) and p^(e) are given by 



(4.4) 



Pn(e) = < 



2 
2 



1- 



1- 



with 



0, 

end, and 



end < e < 



< e < 



(4.5) 



/ + 00 
de/(e)p*(e) 
-00 



2(2/ + 1) 



00 
end I 



|e^d| 



(4.8) 



Her&nJLjUttinger's sum rule or the Fermi-surface sum 
ruleliyO is made use of. 



In Eq. ( [4.4D , c is the weight of the triangular part. 
How ever , it is never a free parameter in this paper. When 
Eq. (2^) is satisfied, our model is nothing but the fiat- 
band model. Such a fiat-band model is simulated by 
p*(e) with c 1/(2/ + 1). When interband transfer terms 
£^^_/(k) are small, on the other hand, all the 2/ + 1 bands 
are degenerate. Such a model is simulated by p*(e) with 
c = 1 or c = 0. 

It follows that 



P,(0)=cV:(0)A„ 



+ (1 - c)v:(o)A 



^wd I 



+c(i-c)p;(o) 



2e, 



nu A I \(^wd\ 
A.1 



2|end| 



An 



I end I 



with 

3 1 — X 



(1 - x) [(1 - x) + 3x] 



and 



Anw{x) = Inx + (x + 1)^ In ( 1 + - ) - X. (4.12) 



(4.9) 



In(l-x)^, (4.10) 



Ay,{x) = -2 [x\nx + (1 - x) ln(l - x)] , (4.11) 



■ 5 ^wd ^ e ^ c-ij 



0, 



e < e 



wd-) ^wu 



< e 



(4.6) 



with Wu, = e^u — ^wd' For the sake of simplicity, our 
examination is restricted to a realistic case, where a fiat 
band lies within broad bands: e^d < < Cnm ^wd < < 

Cyju, and Wyj ^ Wn so that 



p:(o)»/o;(o). 



(4.7) 



It is also assumed that the quasi-particle density of states 
is normalized in such a way that 



+ CO 



1. 



dep*[e) 

) 

Then, the carrier density per site is given by 



FIG. 4. A solid line shows An(x), a dotted line Aw{x)^ 
and a broken line Anw{x). 



The instability condition ( [4. 2D is more easily satis- 
fied for smaller Ps(0). Fig. ^ shows An(x), Aw{x)^ and 
Anu^{x). Note that ^^(0) = ^^(1) = 0; A^(0) = 
^4^(1) = 0; Anw{^) = and Anw{x) increases only loga- 
rithmically with increasing x in such a way that 



(4.13) 



for X ^ 1. When the chemical potential is around one 
of the band-edges, ferromagnetism is favored. On the 
other hand, An{x) and take their largest values 

at X = 1/2; An{l/2) = -(2/3) + (8/3) ln2 = 1.1873 • • • 
and A^{l/2) = 21n2 = 1.38629 • • •. When the chemical 
potential is around the band center, ferromagnetism is 
disfavored. This is one of the reasons why the band-edge 
condition is desired to be satisfied for itinerant-electron 
ferromagnetism. 



B. Flat-band model 

Consider the limit of w^/wn +oo or 

P*JO)/pm - +0 (4.14) 
in the case of nonzero c (c > 0). In this limit, 

P,(0)/p*(0) = cAn {\end\/wn) • (4.15) 

Note that c :^ 1/(2/ + 1) is smaller than unity. This is 
one of the reasons whyjerromagnetism is favored in the 
flat-band model. Eq. (|4.15| ) vanishes in the limit of 



or 



/ +00 



\end\/Wn +0. 



(4.16) 



(4.17) 



Eq. ( [4.17[ ) is nothing but the band-edge condition. Then, 
the instability condition (4^) with AXs = becomes 
quite simple as 



1 - {4>s/4>^) < 0. 



(4.18) 



In the Kondo problem, (ps/^^-f is called the Wilson ra- 
tio. Call it the local Wilson ratio in this paper. In gen- 
eral, it is close to unity in the low density limit of n ^ +0 
and is larger than unity for n > 1. The carrier density is 
simply given by 



n = 2(2/ + 1) (I 

^wd I / ) • 



(4.19) 



When it is close to unity, n :^ 1, charge fluctuations are 
substantially suppressed, so that and ^o, 

as is examined in Appendix Then, it follows from 
Eq. (|0) that 



30s 



4/(/ + l)(^o 



2(2/ + 3) (2/ + l)(2/ + 3)' 

In the case of U = U' and J = 0, it follows that 
so that the local Wilson ratio is as large as 



(4.20) 



1 



5/ + 3 



8/2 + 11/ + 3 



> 1. 



(4.21) 



Here, (j)c = ^ is assumed. When U > or J > 0, 
it follows that > (po- The local Wilson ratio is larger 



tha n Eq . (4.21) so that the instability condition, Eq. (4.2) 
or ( 4.181) , is satisfied in such a case. As long as the Hund 



coupling is nonzero (J > 0), the ferromagnetic instability 
condition is satisfied in the flat-band model with large 
enough / and d. 



C. Band-edge model 

Consider the low density limit of n ^ +0 in the case 
of c = 1 or c = 0. According to the above argument, 
the chemical potential should be at the band bottom. It 
follows that Ps{0) +0 in the limit of n +0. The 
instability condition is also given by Eq. ( [4. 181 ). 

In the single-band model, \imn^-^o(^(j)s/(f)^) — 1 but 
lim^^+o (^5/(^7) < 1 as is shown in Appendix The su- 
perexchange interaction is antiferromagnetic in general, 
and it vanishes in the limit of U/kBTk +00. The 
single-band model with n ^ and U/ksTK +00 is 
close to a quantum critical point for ferromagnetic insta- 
bility. However, the instability condition is not satisfied. 

In multi-band models, on the other hand, it follows 
that lim^^+oo hmn^+o(^s/^7) = 1 for J = while 
lim^^+oo hm^^+o (^s/^7) > 1 for J > 0. The superex- 
change interaction is ferromagnetic for non-zero J > 
and large enough /. As long as the Hund coupling is 
nonzero (J > 0), t herefore, the ferromagnetic instability 
condition Eq. (4^) is generally satisfied in the limit of low 
densities (n +0), large band multiplicity (/ +00), 
and large spatial dimensions ((i ^ +oo).l3 



V. DISCUSSION 

In the single-band model, the superexchange interac- 
tion is antiferromagnetic and it vanishes in the large limit 
of U. The novel exchange interaction is antiferromag- 
netic for almost half filling, for example, according to 
Fig. ||. Therefore, the ground state must be antiferro- 
magnetic or close to a antiferromagnetic critical point in 
the thermodynamic limit of small hole densities discussed 
in Introduction. Ferromagnetic instability cannot occur 
in this limit. Nagaoka and Thouless's ferromagnetism is 
a singular point as a-function of carrier numbers. 

Kanamori's resultti can also be explained within the 
framework of this paper. Both the flat-band and band- 
edge conditions are satisfied in Kanamori's model. Then, 
the novel exchange interaction is ferromagnetic so that 
the ground state must be a paramagnetic state close to 
a ferromagnetic quantum critical point. One can argue 
that ferromagnetic instability can occur in Kanamori's 
model if a weak ferromagnetic superexchange interaction 
is phenomenologically included; the superexchange inter- 
action cannot be ferromagnetic in any single-band model. 

In tlj^|-|flat-band models considered by Mielke and 
Tasaki,B~l3 different orbits are located at different posi- 
tions within a unit cell. In general, the Coulomb inter- 
action between different positions is much weaker than 
the on-site one. Their models correspond to the case 
of > 0, = and J = in Eq. (|2j) within the 
framework of this paper. When and J are small or 
when U' = J, the derivation of the superexchange in- 
teraction in this paper breaks down, as is implied by 
Eq. (3.44), where its denominators vanish. Once a large 



enough but still small U' {U' > J) between different or- 
bits is introduced, the formulation of this paper can be 
applied to such variant models of Mielke and Tasaki's. A 
strong ferromagnetic superexchange interaction exists in 
such models. This argument implies the existence of a 
strong ferromagnetic exchange interaction or an effective 



strong Hund coupling in the flat-band model of Mielke 
and Tasaki's, which corresponds to the ferromagnetic su- 
perexchange interaction considered in this paper. This 
is another reason why ferromagnetism is favored in the 
flat-band models. 

Arguments on the flat-band condition in Sec. |^ are 
restricted to a rather realistic case, where a flat band 
lies within broad bands. If the hybridization or mixing 
between flat and broad bands is strong as it is in our 
model, the treatment of this paper can be extended to 
the case where a flat band lies outside broad bands. Be- 
cause of crystalline fields in cubic lattices, for example, d 
orbit als split into de and dj multiplets. A model where 
only de or dj multiplet is taken into account corresponds 
to / = 1 or / = 1/2, while a model where all d orbitals are 
taken into account corresponds to / = 2. Although the 
superexchange interaction is definitely ferromagnetic for 
large /'s, whether it is ferromagnetic or antiferromagnetic 
depends rather sensitively on / in the region of small Vs 
such as / = 1 and / = 2, as are shown in Figs. H and |[ 
Not only orbitals around the chemical potential but also 
all other orbitals far from the chemical potential should 
be taken into account for a quantitative discussion on 
the magnitude of the superexchange interaction as well 
as its signs. It is interesting to extend the theory of this 
paper to treat such models with non-degenerate d or / 
orbit als a nd models without the symmetry of Eqs. ( |2.4[) 
and ( ^.10| ). Such extensions seem to be indispensable to 
explain quantitatively itinerant-electron ferromagnetism 
in actual metals or alloys. 

It is assumed in this paper that only electrons or holes 
exist. It is also interesting to extend the theory to treat 
ferromagnetism in semimetals, where electrons and holes 
coexist. 

Because the superexchange was originally derived in 
insulating phases, it is controversial whether or not it ex- 
ists in metallic pha ses. The derivation of this paper is 
based on Eq. ( 3.31 ) together with the Ward identity.L3 
It is obvious that the superexchange interaction exists 
even in metallic phases as long as the Mott-Hubbard 
splitting exists. When localized levels are so deep that 

— /i <C —ksTK^ Eq. ( |3.31 ) is a good approximation 
only for l^^l ^ ksTK- In this case, Gutzwiller's nar- 
row quasi-particle band appear at the chemical potential 
between the lower and upper Hubbard bands. The sys- 
tem is a metal because of Gutzwiller's band. However, 
Gutzwiller's band plays no role in the process of the su- 
perexchange interaction. When localized levels are so 
shallow that — ^ —^bTk^ Eq. ( 3.31 ) is a rather 
good approximation even for t he vi cinity of the chemi- 
cal potential; according to Eq. ( ^.3l| ), the mass enhance- 
ment factor is 0^ = 1 — n[2//(2/ + 1)], which is consistent 
with Gutzwiller's result given by Eq. (|L3|) in the large 
limit of / +00. Gutzwiller's band and the lower Hub- 
bard band merge into a band and the upper Hubbard 
band lies above the merged band. In both the cases, the 
Mott-Hubbard splitting exists so that the superexchange 
interaction exists. 

The superexchange interaction in metallic phases can 
also be understood by the following argument. Be- 
cause the typical time scale of quasi-particles, which is 
0{h/kBTB)^ is much longer than that of pair excitations 
of electrons or holes across the Mott-Hubbard gap, which 
is 0{h/U) or 0{Ti/{U' ± J)), all the exchange processes 



can be finished while a pair of quasi-particles meet occa- 
sionally and part again. As long as the Mott-Hubbard 
splitting exists and Eq. ( ^.26[) is satisfied, therefore, the 
superexchange interaction works between quasi-particle 
pairs that occasionally come close to each other. 



The novel exchange interaction considered in Sec. HI C 
is similar to the double and the RKKY exchange inter- 
actions, each of which arises from that of pair excita- 
tions of conduction electrons. The double and the RKKY 
exchange interactions are essentially the same; in many 
cases the double exchange interaction seems to stand for 
the q = component of the RKKY exchange interac- 
tion. However, the novel exchange interaction is totally 
different from the double or RKKY exchange interaction. 
When the RKKY or double exchange interaction is de- 
rived, the existence of localized spins are assumed. The 
double or RKKY exchange interaction is an exchange in- 
teraction working between localized spins in metals. On 
the other hand, the novel exchange interaction works 
between itinerant quasiparticles. The strength of the 
RKKY or double interactions is inversely proportional to 
the width of conduction bands, while the strength of the 
novel exchange interaction is proportional to the band- 
width of quasiparticles. 

Because the local or intrasite part of the RKKY or 
double exchange interaction can play no role, its intrasite 
part should be subtracted. Then, whether the RKKY or 
double exchange interaction is ferromagnetic or antiferro- 
magnetic depends on the position of the chemical poten- 
tial. The band-edge condition also applies to this case. 
When the chemical potential is at the top or bottom re- 
gion of conduction bands, it is ferromagnetic. When the 
chemical potential lies in the center region of conduction 
bands or the nesting of the Fermi surface is significant, 
it is antiferromagnetic. The situation is similar to that 
for the novel exchange interaction. 

The double or RKKY exchange interaction arises from 
the second-order perturbation in the s-d exchange. When 
the s-d exchange interaction is ferromagnetic, the volume 
of the Fermi surface is simply given by the number of con- 
duction electrons. In the band-edge model, the RKKY 
exchange interaction is ferromagnetic and metallic fer- 
romagnetism is possible. Ferromagnetism caused by this 
type of the double or RKKY exchange interaction is noth- 
ing but ferromagnetism examined by Zener.E2l It should 
never be called itinerant-electron ferromagnetism, but it 
is local- moment ferromagnetism in metals. 

When the s-d exchange interaction is antiferromag- 
netic, localized spins must be screened by the Kondo ef- 
fect. The analogy between the periodic antiferromagnetic 
s-d model and the periodic Anderson model implies that 
when higher-order processes in the s-d exchange interac- 
tion are considered the double or RKKY exchange inter- 
action should be replaced by a similar one to the novel 
exchange interaction studied in this paper. The analogy 
also implies that the volume of the Fermi surface must be 
given by the number of total electrons or by the sum of 
those of conduction electrons and localized spins. When 
the volume of electron or hole Fermi surface is small, this 
exchange interaction must be ferromagnetic. For exam- 
ple, assume that there is only a single conduction band. 
Denote the average number of conduction electrons per 
site by Uc- When the number of localized and screened 
spins is taken into account, ric = I corresponds to the 



total electron number 2 per site. This exchange inter- 
action must be ferromagnetic for nc — 1, because the 
band edge condition is satisfied. We speculate that if the 
critical temperature is much lower than the Kondo tem- 
perature ferromagnetism caused by this exchange inter- 
action must be physically the same asJtinerant-electron 
ferromagnetism studied in this paperO 
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VI. SUMMARY 

Itinerant-electron ferromagnetism in multi-band mod- 
els is examined within a theoretical framework of the 
competition between the two leading-order effects inl/d: 
magnetic instability caused by Weiss' mean fields or ex- 
change interactions and the quenching of magnetic mo- 
ments by the Kondo effect. 

Two exchange interactions are responsible for ferro- 
magnetic instability. One is the superexchange interac- 
tion, which arises from the virtual exchange of pair ex- 
citations of electrons or holes across the Mott-Hubbard 
gap. When the Hund coupling is strong enough and the 
band multiplicity is large enough, the superexchange in- 
teraction becomes ferromagnetic. The other is the novel 
exchange interaction arising from the virtual exchange of 
pair excitations of quasi-particles, which is also responsi- 
ble for the Curie- Weiss law. Its strength is proportional 
to the bandwidth of quasi-particles, and its signs depend 
on the shape of the density of quasiparticle states and 
the position of the chemical potential. This exchange 
interaction is ferromagnetic, when the fiat-band condi- 
tion is satisfied so that the density of states has a sharp 
peak around the chemical potential within or outside a 
broad band or broad bands or when the band-edge con- 
dition is satisfied so that the chemical potential is at one 
of the band edges. When the electron or hole density is 
close to unity per unit cell, charge fluctuations are sub- 
stantially suppressed so that the local Wilson ratio is en- 
hanced. Because spin fluctuations are more enhanced 
than orbital fluctuations in the presence of the Hund 
coupling, the local Wilson ratio is also enhanced. The 
enhancement of the local Wilson ratio leads to the en- 
hancement of the novel exchange interaction. When the 
sum of the superexchange and the novel exchange inter- 
actions overcomes the quenching of magnetic moments or 
the suppression of magnetic instability by local and inter- 
site spin fluctuations, a ferromagnetic state is stabilized 
at low enough temperatures. 

In conclusion, the following three conditions should 
be satisfled in order that a ferromagnetic state is sta- 
bilized at low enough temperatures: (i) There are two or 
more than two bands, (ii) the on-site Hubbard repulsion 
is much larger than the energy scale of local quantum 
spin fluctuations, and (iii) the Hund coupling is strong 
enough. Some or all of the following conditions should 
also be satisfled: (iv) There is an almost dispersionless 
band within broad ones, (v) the chemical potential is at 
the top or bottom of bands, and (vi) the electron or hole 
density is close to unity per unit cell. 
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Consider Ha given by Eq. ( p.l3| ). In this Appendix, it 
is assumed that the hybridization energy deflned by 



(Al) 



does not depend on e as well as m and the conduction 
bandwidth 2D is much larger than [/, and J. When 
U = U' = J = 0, the density of states of localized elec- 
trons is given by 



1 



Po(e) = - 



(A2) 



When the Friedel sum rule or the Fermi-liquid relationcil 
is made use of, the density of states at the chemical po- 
tential for non-zero [/, U' and J is simply given by 



(A3) 



as a function of n = Xl^cr(^^cr)- study in this Ap- 
pendix is also restricted to n < 1. 

Introduce inflnitesimally small external flelds given by 
Eq. (3.1) with Hgo = 0, and deflne 



ema = Ca - A/i - ^ags/J^sHg - mgoliBHo. (A4) 

The ground^tate wave function, <I>, is expanded in such 
a way thattll 



$ = 



mkcr m'pa' 



f t 



with |0) the Fermi vacuum, and it satisfles 
[E-Ha- Hext] ^ = 0. 



|0), (A5) 



(A6) 



In Eq. ( |A5D , r[(^^)(^/^/).(p^/)(ka)] is an antisymmetrized 
coefficient deflned in such a way that 



r[(mcr)(m'cr');(pcr')(kcr)] 

— r(^cr)(m'crO;(kcr)(pcr') +r(^/cr/)(mcr);(kcr)(pcrO • (^7) 



r(mcr)(m'cr');(pcr')(kcr) ^ {m' a'){rn(j)\{ji(j'){\<i(j) 



Eqs. (A5) and (|A6| ) give a set of equations for F's: 

- ^ <rr(ma);(ka) = 0, (A8a) 
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"77^ ^m'pr[(ma)(m'aO;(paO(ka)] = 0, (A8b) 
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--y=^{ma);{\c(T) H = 0, (A8c) 
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••• = 0, 



(A8g) 
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m — a ^m' a 



U' 



+ -^r[(^cr)(m'-cr);(p-cr)(kcr)] H = (ASh) 



and so on. In Eqs. (|A8^ , (|A8g|) and (|A8h|) , m ^ m' . 

When s mall irr eleva nt terms are ignored, it follows 
from Eqs. (|A8a|)- ([A8h|) that 



1 
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with 



1 



U + €a- H 
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31 I 
1 
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, (AlO) 



= ea - 5*A/i - -(jglnBHs - mg*iiBHo- (A12) 
Here, (/*, 5* and (/* are renormalized 5 factors: 
<?: = l + ^|^[i?(C/) + 3?i?(C/'-J) + + J)], (A13) 
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(AM) 



and 



^ = 1 . 
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with 



R{x) 



4(2/ + l)A 



-/i D-\-x- 



(A15) 
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Eq. ( |A15D should be discarded for / = 0. Note that > 
9t/9s > 9l/9o > ^'^ 9t/9s = 91/ 9 o when U = U' and 
J = 0, while ^*/^s > 91/ 9o when or J > 0. Even 

in the absence of the spin-orbit interaction, the g factor 
in the spin channel is enhanced in multi-band models. 

When {rima) < 1, F{E) is approximated by F{E) 
— 1/E. The electron number is given by n = —dE/dA/i. 
The charge susceptibility is calculated so that 



x:(o) = 



dn 
dA/1 



{9 c) 



.2 



2(2/ + l)A 



n) {All) 



as a function of n. In similar manners, spin and orbital 
susceptibilities are calculated so that 



ga) 2{2l + l)Agl-n' 



(A18) 



The expansion coefficients appearing in the single-site 
selfenergy are given by 



<t>c = {9c) {9c -^)- 



sin {TTUma) 



for the charge channel, and 



1 {-T^rimaf 



9a J 9i Sm^TTTlma) 



(A19) 



(A20) 



for spin {a = s) and orbital {a = o) channels. Note that 
> 0s ^ 0o > 1 in the low density limit of n ^ +0 
although (j)s > (l)o ^ I ^ 4>c ^OT n 0:^ 1. 
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